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Abstract 



> 

■ The recently proposed strategy for studying the equihbrium thermodynamics of the glass 

^ . phase using a molecular liquid is reviewed and tested in details on the solvable case of the p-spin 

model. We derive the general phase diagram, and confirm the validity of this procedure. We 
point out the efficacy of a system of two weakly coupled copies in order to identify the glass 
transition, and the necessity to study a system with m < 1 copies ('clones') of the original 
\ problem in order to derive the thermodynamic properties of the glass phase. 



e 

^ 1 Introduction 



o 

^ ' It has been proposed recently that the thermodynamic properties of the glass phase can be 

deduced from those of the corresponding liquid phase, computed above the glass transition but 
^ , for a 'replicated' liquid involving m species with m < 1, forming then a molecular hound state 

\ This scheme is appealing in that it provides a general strategy for computing the equilibrium 

properties of the glass phase|§]. This whole strategy is supposed to be valid whenever the 
thermodynamic transition is due to an 'entropy crisis' scenario ||3|, ^, ^, ^ |^, ^ . This basically 
means that the configurational entropy defined as the logarithm of the number of accessible 
metastable states of the system, which is extensive above the transition temperature, becomes 
zero below the transition temperature. 

The archetypes of such systems, and the only ones for which one can get a full control so 
far, are provided by the so-called p-spin models, some infinite range spin glasses involving p > 2 
spin interactions, which have been already much studied. In this paper we want to test the 
general strategy of on these exactly solvable cases. We show explicitely how the properties 
of the glass phase in these systems can be deduced from that of the replicated liquid, and we 
point out some subtleties of the phase diagram of this replicated liquid, which should be kept 
in mind when one deals with the real systems of interacting particles. 

This paper complements the work by Franz and Parisi [^] which performed a similar study 
of the phase diagram of two coupled systems when one of the system is thermalised in the low 
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temperature phase. Here we concentrate instead on the case where we use an arbitrary number 
m of coupled systems which are all identical, as used in Q (The case m = 2 was alluded to in 
PI). 



2 The basic strategy 

2.1 Identifying the glass phase 

Let us assume that there exists a glass phase in some system. How can one distinguish the 
glass phase from the liquid one, in the framework of equilibrium statistical mechanics? A 
simple answer, developed initially in the spin glass context 11, 13], is the following. 
Take two copies, with a small attraction between particles in the two copies, of strength e. 
The attraction is short range (range less or equal to the typical interparticle distance), and it 
couples all particles of system 1 to all particles of system 2. Look at the correlations between 
the two copies when e — > (limit taken after the thermodynamic limit). If there remain some 
correlations in this double limit, this means that the system is in a solid phase (a glass or a 
crystal); otherwise it is in the liquid phase. This general technique is a very robust one which 
is able to identify the glass phase in very diverse systems, ranging from spin glasses to directed 
polymers and to structural glasses [10, 11, |l2|, U, |l|. 



In structural glasses this transition is associated with a spontaneous breaking of a transla- 
tion symmetry. For e = 0, the Hamiltonian is invariant by the relative global translation of 
the particles in the second system relative to the particles in the first one (assuming periodic 
boundary conditions). This symmetry is spontaneously broken in the solid phase. Therefore 
one expects that there is no way in which one can connect continuously the liquid phase to the 
glass phase without encountering a phase transition. We shall check this property explicitely on 
the case of the p-spin mean field spin glass: we shall work out the phase diagram in the plane 
e, T, showing the existence of a glass transition temperature at all e. 



2.2 Nature of the glass transition in the discontinuous case 

This identification of the transition from the study of two coupled system does not allow to 
reach the properties of the low temperature phase from a study of the high temperature one, 
precisely because of the existence of a phase transition. It turns out that, for a specific class 
of systems, there exists a more subtle generalization of this approach which uses an arbitrary 
number m of coupled systems and allows to study the glass phase. These are the systems for 
which the glass transition is of a discontinuous type, called 'one step replica symmetry breaking' 
(one step rsb) in the spin glass language. Basically this situation is encountered whenever there 
exist some low-lying states of the system which are well separated in phase space and have 
independent free energies [^], a situation which has been argued to give a good description 
of the structural glass phase 1^]. More precisely, the extensive part of their free energies are 
equal, but the fluctuating, non-extensive part of their free energies are not correlated from state 
to state. One may expect that such systems will be generic in the sense that this situation 
corresponds to a universality class for extreme event statistics [^]. Proving whether a given 
physical system actually falls in such a class may be quite difficult, but a useful approach, 
adopted in [Q], consists in assuming that it does, computing the properties of the system within 
such an assumption, and comparing them to experimental and numerical data. The success 
of this strategy, and the number of experimental obsevations that it explains and relates, have 
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added confidence to the validity of the old conjecture according to which fragile structural 
glasses [§] fall into this category |^ Q . 

In this paper we shall concentrate on the mean field spin glasses with p-spin interactions 
{p > 2) which have been shown to belong to this class of systems, both for Ising spins and for 
spherical spins 0, |l^, |l^ (Several other sytems have been proven to belong to this class, 
including some spin glasses without quenched disorder |19||). We shall briefly review some of 
their properties which are crucial for understanding the nature of their glass transition and how 
the glass phase can be studied. These systems possess an infinity of metastable states, which we 
shall label by an index a. Each state is characterized [pO| by the value of the local magnetization 
on each site, mf = {si)a- The states, and their free energies, can be defined as solutions of 
TAP equations. A very important quantity is the number M of such metastable states with a 
given free energy density /. It turns out that, in some range of free energies, this number is 
exponentially large in the number of degrees of freedom. One can thus define a corresponding 
entropy, called the configurational entropy in the glass litterature |^ and the complexity in the 
spin glass litterature |21, ^ |6|, defined by 



M{f,T,N)^eMNSc{f,T)), (1) 

The number of metastable states vanishes outside of the region fminiT) < f < fmaxiT), and 
the configurational entropy Sc{f,T) goes to zero at fminiT). 

Let us discuss the properties of such a system at thermal equilibrium. We call the free 
energy density of state a. At low enough temperatures the total free-energy-density of the 
system [fs) can be well approximated by the sum of the contributions to the free energy of each 
particular minimum: 

Z = eM-PNfs) ~ Yl ^M-mfa). (2) 

a 

For large values of we can write 

/fmax 
dfexpi-NiPf-Scif,T)). (3) 
J min 

Using the saddle point method we find that 

fs = ujn<^U) = r-TS,ir,T), (4) 

where 

^f)^f-TS,{f,T). (5) 

This formula is quite similar to the usual formula for the free energy, i.e. / = mmE{E — TS{E)), 
where S{E) is the entropy density as a function of the energy density. The main difference is 
the fact that the total entropy of the system has been decomposed into the contribution due 
to small fluctuations around a given configuration (this piece has been included into /), and 
the contribution due to the existence of a large number of locally stable configurations, the 
configurational entropy. 

Calling /* the value of / which minimizes *!*(/), we have two possibilities: 

• The minimum of <I> lies inside the interval and it can be found as the solution of the 
equation (3 = dSc/df. In this case we have 

<i> = r-Ts:, s: = s,{r,T). (6) 
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Figure 1: Qualitative shape of the configurational entropy versus free energy. The whole curve de- 
pends on the temperature. The saddle point which dominates the partition function, for m correlated 
clones, is the point /* such that the slope of the curve equals m/T (for the usual, uncloned, system, 
m = 1). If the temperature is small enough the saddle point sticks to the minimum / = fmin and 
the sytem is in its glass phase. 



The system may stay in one of the many possible minima. The number of accessible 
minima is exp(A^S'*) . The entropy of the system is thus the sum of the entropy of a 
typical minimum and of S*, which is the contribution to the entropy coming from the 
exponentially large number of metastable configurations. 

• The minimum of <I> is at the extreme value of the range of variability of /: it sticks at 
/* = fmin and the total free energy is $ = /mm- In this case the contribution of the 
configurational entropy to the free energy is zero. The different states which contribute to 
the free energy have a difference in free energy density which is of order (a difference 
in total free energy of order 1), and these free energy differences are uncorrelated and 
exponentially distributed ||l5| , |20[| , corresponding to Gumbel's distribution of extremes 
@- 

Therefore the glass transition of this system is a thermodynamic transition associated to what 
is sometimes called an 'entropy crisis': the entropy associated with the number of metastable 
states which contribute to the partition function vanishes below the transition temperature. The 
transition is of second order from the thermodynamic point of view, having a continuous entropy 
and a jump in specific heat, but it is associated with a discontinuity of the order parameter 
We call such a transition a 'discontinuous glass transition'. 



2.3 Cloning 

Clearly the equilibrium thermodynamics of this class of glasses is fully described by the knowl- 
edge of the configurational entropy as a function of temperature and free-energy-density. Monas- 
son 1 23 1 proposed the following method to compute this function. 1 Let us consider m copies of 
the system which are constrained to stay in the same minimum. We shall discuss below how 
one can achieve this constraint, but let us first discuss the physics of this constrained system. 
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Figure 2: General shape of the free energy per particle 0(m, T) of the cloned system, versus the 
number of clones m, at temperatures below the glass transition. The critical value m = mc{T) is 
identified as the first point where the slope d(j)/dm vanishes when one increases m. 



Its partition function is basically the Laplace transform of the number of states: 



Z^= /^"d/e-^t-^^-^^^^'^)] . (7) 

^ fm 



The free energy per spin, 



(Pirn, T) = --i- log Zm « min( / - -S,{f, T)) , (8) 
pmJS f m 

allows to compute the configurational entropy Sc{f,T). Because of the concavity of the Sc 
versus / curve, its largest possible slope is 

«m^^(/™n,T) . (9) 

Performing the minimisation in (P), we find that the states which contribute to the partition 
function have a free energy density f*{m,T) given either by 

^^(/•,T) = = (10) 

if m is small enough (so that the slope at the origin, a{T), be larger than m/T), or otherwise 
they are just the low lying states with /* = fmm (see fig. ||). Therefore, varying m at a fixed 
temperature T < Tc, we have a phase transition at a critical value of m = mc(T) which is 
smaller than one. 

The general shape of the free energy (f){m,T) is the one shown in fig.|2[ For m > rriciT) the 
free energy is constant (p{m,T) = fmin{T)- For m < mc(T) this function allows to compute the 
configurational entropy through the implicit equation: 

d[m(j){m, T)] 



dm 
m? d(l){m, T) 
T dm 



f 

Sc ■ (11) 
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Increasing m, the transition point rric is determined by the vanishing of the slope ^"^g™'"^^ • 

The following two remarks, which are trivial consequences of the simple nature of the phase 
transition in these discontinuous sytems, are crucial: 

• For m < mc{T), the system is in its high temperature phase, even if T < Tc. The free 
energy (j){m, T) is that of the high temperature phase of a system of m constrained copies, 
a liquid made of molecules built up from m atoms. It can thus be evaluated by any of 
the methods giving a good description of the high temperature phase, including series 
expansion, simulations etc... 

• The transition at m = mc{T) is such that the free energy is constant for m > mc{T) (this 
is a consequence of the transition being due to the sticking of the saddle point in / at the 
lowest value fmin)- Therefore the free energy in the glass phase at a temperature T < Tc 
is given by F{T) = (j){m = 1, T) = <j){mc{T),T). But the free energy is continuous through 
the transition line, and therefore (j){mc{T),T) is the free energy of the high temperature, 
'liquid', phase. 

This is the method which has been used in Q to compute the thermodynamics of the glass. 
It involves dealing with a cloned system made up of m identical clones constrained to be in the 
same state, and to do this for any m, in such a way as to be able to analytically continue the 
results to m < 1. This continuation bears some similarities to, but is different from, the one used 
in the replica method in order to compute quenched averages. As we shall also need to use these 
more standard replicas, we shall call the m copies introduced here and in |2^, |l|] some 'clones' of 
the original problem (we keep the word copies for the case when m is an integer larger or equal 
to one). One can introduce the clones even for systems which have no quenched in disorder, 
like structural glasses. The reader will have noticed that the above discussion is a general one 
which does not rely on the details of the p-spin mean field models; in fact the procedure of 
computing the free energy of the cloned and constrained system, (j){m), in the liquid phase, can 
be used also in some non mean field systems [Q, in which it is more difficult to come up with 
a definition of the metastable states, since these involve fixing some state-dependent boundary 
conditions in a self consistent way. The replicas introduced in disordered systems, which are 
used to handle the average over quenched disorder, can be introduced on top of the clones if 
needed: in the case of systems with quenched disorder, clones can be replicated! As we shall 
see, the computation with the fixed number m of coupled clones in the high temperature phase 
is equivalent for disordered systems to a 'one step rsb' computation. This has been first found 
in some computations of neural networks |27] and is expected on general grounds [p^j . 

We would like to also comment on the alternative method which has been proposed inde- 
pendently for the study of the glass phase ||2^. It also introduces some clones of the system, 
but it uses one reference system which is supposed to be thermalised in the glass phase, and 
then studies the free energy of some other systems, constrained to be at a fixed distance of 
the reference one. This is a powerful method which has the advantage of allowing to study 
of metastable states in the high temperature phase, as has been demonstrated in ||2^, It 
has the drawback of requiring to have one system thermalised in the glass phase, which is not 
easy to achieve. The strategy of cloning discussed here, which has been used in |^, || has the 
advantage of involving only computations in the liquid phase, but its scope is probably limited 
to the study of equilibrium thermodynamics. 
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2.4 Two coupled copies 

We have seen in the previous section how a set of m clones constrained to stay in the same 
state can allow for the computation of thermodynamical properties of glasses. The question we 
address now is how to achieve this constraint. In order to try and constrain the m clones to 
be in the same state, we shall introduce an attractive coupling between the various clones. For 
spin systems, calling the spin at site i in clone a, this coupling could be for instance: 

TV 

E E<4- (12) 

\<a<h<m i=l 

The detailed structure of this coupling is to some extent irrelevant, but it should satisfy the 
following constraints: it is extensive, in that it contributes to the energy at order N, and it is pro- 
portional to e. What we shall be mostly interested in is taking the double limit limf^ — limTv — ^co. 
But it will also turn out to be useful to study the phase diagram of the coupled system for 
general e. This we shall do in the next section for the p-spin models, but let us first try to get 
a qualitative understanding of the phase diagram in general, based only on the existence of the 
configurational entropy. The various parameters in our problem are the number of clones m, 
the temperature T, and the coupling between clones e. 

Let us first consider the case of two coupled copies, m = 2, in the limit of small e. There 
can be several contributions to the partition function: 

- When the two systems fall into the same state, we have: 

df eN[Sc{f,T)-2f/T+eqi/T] ^ ^-2N4>=/T (^^3^ 

:n(T) 

where qi is the typical self overlap of a state. 

- When they are in different states, we have: 

Z^c^ [ df e^[25c(/'^)-2//^l = e-^^-^^/^ (14) 

fmin {^) 

There exist several temperature regimes: 

- If T < Tc(0), both contributions are dominated by the boundary / = fmin, giving (f)= = 
fmin — and 4>jt = fmin- Therefore for any positive e, we have (p= < (p^ and the leading 
contribution is due to systems falling into the same state. The system is in a correlated phase 
because the overlap between the two copies is non zero, and it is in the glass phase because the 
free energy sticks at fmin- 

- When Tc(0) < T < T2, where T2 is the temperature at which a{T2) = 2/T2, we have 
a{T) < 2/T. Thus the contribution Z= is dominated by the boundary and ^= = fmin — e9i/2. 
The contribution is dominated by a non trivial sadle point /* such that: {dSc/df){f*) = 
1/T. The corresponding free energy is (j)^ = f* — TSdf*)- Depending on the value of e, the 
minimal free energy will be either (f)=, in which case the system is in a correlated glass phase, 
or (f)^, in which case the system is in an uncorrected liquid phase. Because of the concavity 
of the Sc{f) curve, /* — TSdf*) < fmin and therefore the 'liquid' saddle point dominates at 
e small enough. The critical value of e locating the transition to the glass phase, with the two 
systems sticking together, is 

ec = ifmin-f* + TSdf*))/qi (15) 
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This whole discussion holds only at small e, and in this region one can locate the behaviour 
of the transition curve in the neighborhood of the point T = Tc, e = 0, by using the following 
expansion of the configurational entropy at / > fmin- 
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We have then: 



Scif) = a{T){f - Un{T)) - ^(/ - Un{T)f + 0({f- Un{T)f) . (16) 



/* ^ + ^ {a{T) - l/T) (17) 

In the neighborhood of the point T = Tc,€ = 0, the generic behaviour is that a{T) — l/T is 
linear in T - rc(0). Writing a{T) - l/T ~ (/(T - Tc(0)), we get finally: 



This discussion suggest that, in the limit of small coupling (e — > 0), the system of two coupled 
copies will have a transition to the glass state at a temperature which depends continuously 
on e, with a square root singularity. This transition is clearly identified by a first order jump 
in the correlation between the two systems, defined here as the overlap. Such a system with 
two copies may thus be a good candidate for the numerical determination of Tc{0), as already 



discussed in [25|. 



Let us also notice that there might in principle exist another phase of the system, with the 
two copies correlated, but being in the liquid phase. (This is the phase of the molecular liquid). 
This can exist whenever T2 < T: The contribution is then dominated by a non trivial sadlle 

point such that: (dSc/df) (^f^ = l/T. This is also the case of the contribution Z=, with 
a saddle point fl such that: {dSc/df){fl) = 2/T. Concavity implies that (p^ < at small 
enough e. Therefore the system is in its uncorrelated liquid phase for e small enough, at all 
temperatures T > T^. in a system with two coupled copies one does not see the correlated liquid 
phase at small e. If one could use these small e computations to induce some conclusions for 
finite e, they would hint at a transition, in the liquid phase, between the uncorrelated one at 
small e and the correlated one at larger e. 

In fact for very large e the two coupled systems will be in the same configuration, giving a 
partition function Z = J exp{—2PH). Therefore this system will have a usual glass transition 
at a static temperature: 

lim r,(e) = 2r,(0) (19) 

e — ^00 

and this transition will be a discontinuous glass transition, therefore it will have no entropy 
jump. 

The explicit computations in the p-spin models will confirm the various aspects of the phase 
diagram of the two copies system that we have predicted here. We expect that these are general 
features which will be found in all the glass phases. 



2.5 Coupled clones 

Let us now study the main aspects of the phase diagram that one can expect when there is a 
general number, m, of clones. As we have seen above, the system can be, for small e, in any of 
the following three possible phases: 

• A correlated glass phase in which the free energy is (p = fmin — {m — l)€qi/2. 
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• A correlated liquid phase in which the free energy is the saddle point of = / — 
(T/'m)Sc{f,T) — (m — l)egi/2, and the extremum value /* is greater than fmin- 

• An uncorrelated liquid phase in which the free energy is the saddle point of = / — 
TScif, ^) ^^'^ the extremum value /* is greater than fmin 

One has to distinguish two regions in m space. 

If m > 1, the discussion is very similar to the one we did above for the m = 2 phase, and 
we expect that at small e the system will have a transition from an uncorrelated liquid phase 
at T > Tc to a correlated glass phase at T < Tc. 

If m < 1 the situation is different, basically because the interesting part of the phase diagram 
now lies in the T < Tc region. When T < Tc, the uncorrelated liquid cannot exist because 
a(T) > l/T. Therefore one expects a transition from the correlated liquid to the correlated 
glass, at a temperature such that a{T) = m/T. This is precisely the transition which was 
discussed in the previous section. So our conclusion is that the e coupling, for small e and 
T < Tc, is able to polarize the system in the same state provided m < 1, and thus it realizes 
the constrained system that was discussed in sect. and allows for the computation of the 
glass thermodynamics from liquid computations. 



3 Explicit computations with the ^-spin model 

We now move on to some explicit computations of the phase diagram of the cloned p-spin model. 

3.1 The usual replica solution 

We work with a spherical spin glass model involving N spins Si with the spherical constraint 
J2i Sf = N and an energy given by: 

H = - Jn■■■^,Sn.■■S^, (20) 

ii<...<ip 

where the random couplings Ji^...ip are independent gaussian random variables with zero mean 
and variance p\/{2NP~^). In the explicit computations below, and for all the numerical results, 
we shall use the value p = 4; however all the systems with p > 2 have the same qualitative phase 
diagram. We shall not reproduce here the various techniques which have already been used to 
study this system, but just state some results of interest for our discussion, referring the reader 
to the original litterature [0, The replica solution shows the existence of a transition 

from a replica symmetric phase at T > Tc to a phase with one step rsb. The breaking involves 
a parametrization of the overlap matrix between n replicas, Qab = (1/-^) which has 

a one step rsb Parisi structure ||2^ with a diagonal element equal to one, the elements inside 
diagonal blocks of size x equal to qi, and all the other elements equal to zero. The free energy 
is 

F = -^[l-{l- + 1=^ log(l -Qi)-^ log[l - (1 - x)q^] (21) 

and should be maximised over qi and over x in the interval [0, 1]. 

The high temperature phase is obtained when x ^ 1. In the high temperature phase qi = 
and the free energy is just —(3/2. 

The parameter qi is interpreted as the Edwards-Anderson (EA) order parameter, charac- 
teristic of the degree of spontaneous ordering inside each state: qi = (1/-/V) X]j(5'j)a- The low 
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Figure 3: Replica solution of the p=4 model at equilibrium. We plot the Edward-Anderson order 
parameter qi (dashed line) and the one step rsb parameter x (full line) versus the temperature. The 
transition is at the temperature T = .503 



temperature phase has an EA order parameter qi ^ 0, and qi jumps discontinuously at Tc, 
although the transition is thermodynamically of second order, with a continuous entropy and 
energy (this is possible because x — > 1 at the transition | ^]). The parameter x is characteristic 
of the distribution of free energies of the low lying states |15|, pC|] . 

The result for p = 4 is summarised in fig.^, which plots the value of the breakpoint parameter 
X and the EA parameter versus temperature. The transition temperature, identified from x = 1, 
is Tc = .503. The solution with x > 1 is interesting for the study of metastable states and barriers 
[18 1, but here we stay within the framework of the equilibrium measure where it does not play 
any role. 



3.2 Two coupled copies: exact solution and phase diagram 

We now turn to a coupled system with m = 2 identical copies of the original problem, two sets 
of spins Si and 5^', characterized by the energy: 

H = - E Jii-ip + S'i-r--S'i^ - e^SiS[ (22) 

ii<...<ip i 

The thermodynamics of a cloned p-spin model with some coupling between the clones has already 
been studied in |Q in order to understand the structure of metastable states and barriers. We 
shall borrow from this paper the form of the free energy of the coupled system, computed with 
the replica method (replicated clones), within the one step rsb formalism. The replica formalism 
introduces n replicas of each of the two spin systems. The free energy is expressed in terms 
of a self-overlap matrix between n different replicas of the spins Si , Qab = (1/-^) J^ii^i^i) = 
(1/A^) J^iiSi^ S'^) . It also involves the correlations between the replicas of the two sets of spins, 
Pab = {^/N)J2i{'SiS'^). Within the one step rsb solution, which ws shown to be correct in 
[]l8| , the matrix Q is parametrized exactly as in the previous case by the value of its elements, 
qi, inside the diagonal blocks of size m, while the matrix P is parametrized by the value of its 
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elements, pi, inside the diagonal blocks of size m, and the value p of its diagonal elements (all 
other elements being 0). In terms of these order parameters, the free energy is given by: 

F= - /?[/(!) + /(p)-(l-x)(/((7i) + /(pi))] 

- ^-[log(l - P + (1 -x)(pi -gi)) +log(l + p- (1 -x)(pi + gi))] 
1 1 — X 

+ ^—^[^og{l-p + pi-qi)+log{l+p-{pi+qi))]-ep (23) 

Where /(r) = r^/2. We have two real copies, so the free energy must be minimized with respect 
to p. But we have n ^ replicas to handle the quenched average, so the free energy must be 
maximized with respect to pi,qi,x [pO|] . 

The interpretation of the various order parameters is as follows. The system may freeze into 
a correlated glass phase. This will be signalled by the Edwards- Anderson-like order parameters: 

91 = ^ T.iS^)a{Sl)a ^ Pi = ^ J2iS^)a{S^)a (24) 

i i 

On the other hand, for any finite e, and both in the high and low temperature phase, the degree 
of correlations between the two systems is measured by: 

P = ^J:iS^Sl) (25) 

i 

The extremisation of the free energy at p = 4 and finite e, performed by a simple grid search 
followed by a minimisation procedure, gives the phase diagram plotted in fig. ^. 

Above the full and dashed-dotted lines, the system is in a liquid phase, where the EA order 
parameters are qi = Pi = 0. The correlation p is nonzero, but it goes to zero in the limit 
of e — > 0. Below the full and the dashed-dotted lines the system is in the glass phase, with 
9i = Pi > 0. The correlation p is nonzero, and it remains finite when e — > 0. Therefore the 
e ^ transition is between an uncorrelated liquid phase at high temperature and a correlated 
glass phase at low temperatures. It occurs at the temperature T'c(O) = .503, with some small 
e corrections which are of order ^/e, as anticipated in our general discussion of sect. This 
glass transition is a first order transition on the portion with small enough e (full line). On 
this portion the various order parameters qi,pi,p,x jump discontinuously at the transition. In 
particular the correlation parameter p is discontinuous on this transition, which should make 
it relatively easy to detect numerically (letting aside metastability effects), and there is also a 
latent heat. The glass transition at larger e (dashed-dotted line) is of the second order type 
from the thermodynamic point of view. It has a jump in the EA order parameter qi = pi, 
but the correlation parameter p and the entropy are continuous at the transition; it is thus a 
transition which is exactly of the 'discontinuous spin glass type', as it is known to occur both 
for the single system e = 0, as well as for the infinitely coupled system e — > oo. 

The first order transition line goes on into the liquid phase, where there is a portion (dotted 
line) which is a first order transition line. This transition has nothing to do with the glass 
transition: on both sides of the line the EA order parameter is qi = pi = 0. It is a transition 
in the liquid phase between a regime of strong correlations between the two copies (large p) 
when e is larger than the transition value, and a regime of smaller correlations, we shall call 
it the correlation transition. Along this line the entropy is discontinuous. However these two 
transition lines (full line and dotted line) are not the analytic continuation one from the other: 



11 



1.00 
0.90 
0.80 
0.70 
0.60 



'0.0 2.0 4.0 

Sqrt(epsilon) 

Figure 4: Phase diagram in the plane T, for a system of two coupled copies of a p-spin model with 
coupling e, and p = 4. Full and dashed-dotted lines: the glass transition temperature Tc, identified 
by a jump in the Edwards- Anderson parameter, which is zero in the high temperature phase, and 
jumps discontinuously to a positive value in the low temperature phase. The full line is a first order 
transition with a latent heat and a jump in the correlation (see text), while the dashed-dotted line 
is a 'discontinuous spin glass' transition with no latent heat, no jump in the correlation, but a jump 
in the EA order parameters. The dotted line is a first order transition line, ending at a second order 
critical point, separating two regions of the liquid with different correlations between the two copies. 




this is clearly seen from fig. |5| which plots the jump of entropy along the two discontinuous 
transition lines. The latent heat is continuous across the triple point when one follows the 
discontinuous transition lines, but its derivative is discontinuous across the triple point. The first 
order correlation transition ends up in the liquid phase at a critical point where the correlation 
transition becomes of second order. This correlation transition is very similar to the familiar 
liquid-gas transition. 

It is interesting to notice how the glass transition is transmuted by going from a pure system 
to a system with two copies. In the pure system it is a 'discontinuous spin glass' transition 
which is difficult to identify, being only signalled by a jump of the EA order parameter. When 
we switch on a small coupling e between the copies, we get a real first order transition with 
an easily measurable jump in p and a specific heat. This was already noticed by Franz and 
Parisi in their study of this same problem, as well as in their study of two constrained systems, 
one being thermalised in the glass phase |25|. One outcome of the present computation is that 
the glass transition is present at all e: there is no way one can go continuously from the high 
temperature phase to the glass phase without crossing a phase transition. This transition can 
be either of the 'discontinuous spin glass' nature when the coupling e is either zero or large, 
or it may be first order when the coupling is small enough and non zero, because it is then 
accompanied by a correlation transition. 



3.3 Studying the glass through the cloned liquid 

We now move to the case in which we have m coupled clones, and study specifically the case 
m < 1. The general solution could be studied from the formulas derived in [I8|, but here 
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Figure 5: Discontinuity of the entropy along the two hues of first order transition (full line and 
dotted line) in fig. ||. The transition becomes of second order (continuous entropy) at the two end 
points: the end point at e = which is the usual glass transition, and the critical point for the 
decorrelation transition in the liquid phase. The entropy jump is continuous along the triple point, 
but its derivative is discontinuous. 



we want to concentrate on the test of the general method described above, which states that 
the thermodynamics of the glass can be studied from a high temperature computation of a 
correlated system. We thus consider the m systems with a fixed correlation between any of the 
m{m — l)/2 pairs: 

1 E = 9 (26) 



and compute the partition function 



^ = / n (n dsfs (Y^istf - iv) ) n ^ E s^s^ - 

exp[- E Jn...,f: Sl...sA (27) 

y il<...<ip a=l j 

As we perform only a high temperature computation, we do not need to introduce the replicas 
in order to perform the quenched average. The simplest method is to notice that the annealed 
average is exact in the high temperature phase (or we can also expand Z in a high temperature 
series). We find a free energy (p = — log Z/ {(iNm) equal to: 

3 1 — m 1 

= - (1 - m)q^] + log(l - 9) - ^ log[l - (1 - m)q] (28) 

This is exactly the free energy of a single, uncloned system, when it is studied through the 



introduction of replicas in (21). Our number of clones m is equivalent to the breakpoint x in 
the one step rsb Ansatz, and the correlation q is equivalent to the matrix element qi in the 
replica method. As explained in the previous section, this could be expected on general grounds 



23 1 . It also implies that the proposed method of the high temperature expansion of a cloned 
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Figure 6: Phase diagram of the cloned p-spin glass with p=4 at small e in the plane m (number 
of clones) - temperature. The diagram is obtained from maximising (if m < 1) or minimizing the 
free energy of m clones, with the high temperature-liquid free energy involving only the correlation 
q between the clones. Above the full line the system is in the liquid phase, below it is a glass. The 
free energy along the transition line at m < 1 can be obtained from a high temperature computation 
in a correlated liquid. As the free energy in the glass phase is independent of m, this allows to 
compute the free energy in the glass phase. The dotted line is a first order transition line between a 
liquid which is coupled, with a finite value of g, below the dotted line, and one where the clones are 
uncoupled, (g), above the line 




system works. In fact if one studies, for m < 1, the free energy (/> using the general method 
outlined in sect. 2^, one gets the phase diagram in the plane m,T shown in fig. ^. The value 
of the critical temperature Tc{m) is the inverse function of the x{T) curved obtained with the 
replica solution, and the thermodynamics deduced from sect. 2^ is automatically correct. 

The interesting point of this study is the structure of the phase diagram of fig. ^ It confirms 
our expectations from the above general discussion: the transition is rather different depending 
on whether m is larger or smaller than one. When m < 1 the glass transition is between a 
correlated liquid and a glass. Therefore the high temperature phase that must be studied in 
order to compute the properties of the glass is that of the correlated liquid . When m > 1 
the transition is between an uncorrelated liquid and a glass. This explains some of the previous 
failures and successes of recent works on glasses using this donation idea. In the transition 
in soft or hard sphere systems was studied as one from an uncorrelated liquid to a glass. This 
allowed to get a good result for the transition temperature, but it did not give the correct 
glass thermodynamics, because this glass thermodynamics should be studied from a study of 
a correlated liquid, not from an uncorrelated one. The correct procedure was performed in 
through a study of a molecular liquid. 



4 Conclusion 

We have tested the general method of |Q] on the solvable case of the p-spin model. It works 
well in the sense that the study of a 'molecular liquid' with m < 1 clones allows to obtain the 
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exact result for the thermodynamics of the glass phase. Obviously the computations which are 
performed with the clones in their high temperature phase are equivalent to the direct 'one 
step replica symmetry breaking' solution of the problem: the clones are basically equivalent to 
replicas . However one should not forget that the clones give us a method to compute the 

thermodynamic properties of glasses even in systems without any quenched-in disorder, where 
there is no way to introduce the replicas. In fact we believe that the method of replica symmetry 
breaking, although it was historically introduced in order to compute some quenched averages, 
involving the logarithm of a partition function, is much more general. At the one step level it 
is able to address any kind of problem (with or without disorder) where the low lying minima 
are uncorrelated, with a distribution which is unbounded but decreases faster than any power 
law |16|. This is very interesting since these problems build up one natural universality class 
for extreme events statistics. The question of whether a given problem (e.g. a given type of 
structural glass) falls into this universality class is a difficult one. If it does not, it means that 
there will be another phase transition which forbids the direct study of the glass from that of 
the cloned liquid (a system with two step breaking could be studied by cloning the clones...). 
This is very probably what happens in the Edwards Anderson model |31]. The case of structural 
glasses, if it is confirmed that it belongs to this universality class of systems which are solved 
by cloning, (the favourable indications are reviewed in |j|]), is thus a much simpler one. Yet it 
would also be interesting to try to carry out a similar method on finite dimensional spin glasses, 
in order to see the predictions e.g. for the ground state energy. 
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